Origin of Fractal Branching Complexity in the Lung 

STEPHEN H. BENNETT 1 , MARLOWE W. ELDRIDGE 2 , CARLOS E. PUENTE 3 , RUDOLF H. RIEDI 4 , 
THOMAS R. NELSON 5 , BOYD W. GOETZMAN 1 , JAY M. MILSTEIN 1 , SHI AM S. SINGHAL 6 , 
KEITH HORSFIELD 7 , MICHAEL J. WOLDENBERG 8 


Correspondence Address 

Stephen H. Bennett 
Division of Neonatology TB193 
Department of Pediatrics 
University of California, Davis, 95616 


Phone: 530-752-9146, pager 916-762-5353 
Email: shbennett@ucdavis.edu 

Data Sources/ Analysislsite: http://neonatology.ucdavis.edu/fractal 


>4J 


J 

/ 


A 


Affiliations-( present and former*) 




1 . Department of Pediatrics 

Division of Neonatofogy 

University of California, Davis 95616 


2. National Institute of Health 

Center for Environmental Health and Lung Biology 

University of Chapel Hill 

104 Mason Farm Road CB#7130 

Chapel Hill, North Carolina 27599 

3. Institute of Theoretical Dynamics and 
Department of Land Air and Water Resources 
University of California, Davis 95616 

4. Department of Electrical Engineering 
Rice University MS 366, 

6100 Main Street 
Houston TX 77005 - .1892 


5. Department of Radiology, 0610 
Division of Physics 

University of California, San Diego 
Medical Teaching Facility, Rm-1 13 
9500 Gilman Drive 
La Jolla, CA 92093-0610 

6. *Laxmipat Singhania Institute of Cardiology 

GSVM Medical College Kanpur 
Uttar Pradesh, India 


7.* Midhurst Research Institute 
Midhurst, West Sussex 
England 


8. Department, of Geography 
Wilkeson Quad 
State Univ. of New York 
Buffalo, NY 14261 USA 


Origin of Fractal Branching Complexity in the Lung 
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Shiam S. Singhal, Keith Horsfield, Michael J. Woldenberg 

Fractal and non-fractal theories predict universal designs for biological branching transport networks whereby energy 
dissipation is minimized. However, they are inadequate to account for adaptation via designs that appear sub- 
optimal. Therefore we have evaluated another theory in the developing lung which is based on bifurcation design and 
the principle of self-organized criticality. The latter is a universal physical theory whereby fractal and multifractal 
patterns of branching complexity self-organize into meta-stable states poised between order and disorder. We find 
bifurcation design in the lung manifests an adaptable fractal dimension, with spatial scaling properties implicating 
structural origins of airway and vascular disease susceptibility. 

Introduction 

Branching transport networks in biology form patterns adhering to design principles that 
are energetically favorable, a condition considered important for species survival(l). Such 
networks are constructed of vessels via a bifurcating structure, where the parent and daughter 
vessels, of diameters d 0 , d , , and d 2 , respectively, satisfy a bifurcation design relationship 
cIq = d\ + dl exhibiting the property of self-organization (2). The bifurcation exponent, x, is said 

to specify a design because its value reflects the fluid dynamic efficiency and energy dissipation 
during transport(3). At the network level, an index of design, X, reflects a power-law relationship 
between vessel number and diameter N d~ x , where x ~ X(A). Although it is not clear to what 
extent branching complexity implicit in a power law contributes to form and function beyond the 
vessel level (4-6), both fractal and non-fractal theories of biological network organization predict 
universal designs with energy dissipation properties at a global minimum or a relative minimum 
(1,2,7). Not accounted for are biological exceptions, such as the fetal and neonatal pulmonary 
arterial circulation where design at birth begins with an unusually unstable form of organization 
for fluid transport(8). A premise of universal theories is that the end design is static, being 
independent of influences from the environment and unaffected by the functional requirements 
during development. However, not fully anticipated in current theory, is that the embryonic form 
of design may metamorphose into a final form during growth and development, in which the early 
form may satisfy quite different criteria than the final, juvenile or adult state(9) and where 
patterns and timing of adaptation during successive developmental stages may be just as 
important for species survival, as the final design(lO). In this regard, present theories of biological 
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design are limited in extent because they do no offer explanations for adaptation or mal- 
adaptation via transitional forms of organization, only the achievement of an idealized final 
teleological design(l,2,7,9). 

In his tome On Growth and Form , D’Arcy Wentworth Thompson envisioned a universal 
principle governing the development and adaptation of form (9) that may have a foundation in a 
recently proposed physical theory called self-organized criticality! 11). This theory postulates a 
universal process of fractal pattern formation for self-organizing fluid transport networks in 
nature predicting adaptable patterns of complexity not accounted for in present fractal and non- 
fractal theories of biological network design(12-13). For example, present fractal theory of 
biological branching form is based upon the principle of self- similarity where the condition of 
x=X, results in a space-filling fractal structure, such that x=X=D is predicted to be fixed, 
emerging early in development^). However, biological systems do not necessarily exhibit a 
singular power-law with a design that is universal (4), nor are power-laws of biological branching 
networks scale-invariant, a condition necessary for self-similarity! 14). At present, the fractal 
architecture of biological networks is an unresolved controversial conjecture, making it unclear as 
to whether any form of fractal dimension characterizes biological network design(15, 16). 
Alternatively, a non-fractal theory predicts that biological networks are unorganized stochastic 
network of bifurcations whose overall function is unaffected by network scaling, being 
determined by geometric properties residing at the vessel level only(l, 2). Design universality is 
characterized by a prevailing network average of x coinciding with a global minimum or relative 
minimum of energy dissipation(3). However, one universal feature of biological networks not 
accounted for a universal theory is the manifestation of a broad non-normal distribution of 
bifurcation designs for which a simple average does not necessarily characterize the population 
uniquely (4,5,8). As large proportions about the average deviate from ideal principles it is not 
known if the phenomenon is simply biological variation or is related to the optimization of some 
alternative global cost function related to bifurcation design (3). Alternatively, aberrations not 
accounted for by universal theories are predicted to be dynamic elements of complexity in a self- 
organized critical state(12, 17-19), In attempting to build upon Thompson’s vision, more vital is 
the notion that self-organized criticality does not end with an ideal teleological design, but begins 
with a universal process underpinning design formation, development and adaptation that would 
lead to possible explanations of diversity and uniformity in biological form. 
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The case for a fractal self-organized criticality model in biological transport networks is 
compelling. First, fractal forms of branching are widespread in nature and are advantageous for 
transport because such forms optimize energy dissipation relative to other patterns(13). Second, 
self-organized criticality ( 11-12) represents a plausible theory of fractal pattern formation of 
import in explaining the process of biological network adaptation(9). Rather than an a priori 
assumption regarding universal fractal design, a self-organized critical state postulates a universal 
organizational configuration of fractal branching complexity residing in a marginally stable state, 
possibly far removed from a global minimum of energy dissipation^ 2). This theory incorporates 
additional patterns of branching complexity, such as multifractal branching (13) and self-affine 
fractal scaling(12,14) where the effect of complexity is to break symmetry and impart network 
scaling properties that spatially modulate forces in self-organization (13). For complexity to 
emerge, criticality requires the presence of a coupled feedback process to ensure a steady state for 
which a system is marginally stable against a disturbance (12). Universal patterns are postulated 
to develop when two types of feedback processes are coupled, one threshold-dependent, the other 
threshold-independent, (17-19). The coupled process is general so that different transport systems, 
biological or non-living, can form universal patterns of complexity by different mechanisms(20). 
Consequently, where non-living transport systems develop over eons of time via a punctuated 
distribution of forces(13), biological pattern development depends upon the critical timing of 
information expression, constrained forces and indirect environmental circumstances (21). Hence, 
for biological networks, self-organized criticality conceives design as a complex adaptive process 
containing a continuum of physical transitional states with universal patterns of complexity 
appearing at all stages of development(22). Consequently, this theory would be expected to 
exhibit locally emergent fractal scaling behavior, giving insight into a network organizational 
process modulating physiological adaptation and mal-adaptation, rather than a universal fractal 
dimension. 

The mammalian lung is an example of a self-organizing biological branching transport 
network where conclusive evidence discriminating between a simple-fractal, a complex-fractal or 
a non-fractal theory is lacking. We will formalize fractal branching complexity, generalizing self- 
similar and self-affine scaling in power-law behavior, along with patterns of branching for fractal 
and non-fractal design in biological trees, using a dimension defined according to bifurcation 
design, x. From this formulation, alternative theories of biological network organization will be 
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tested in the transport vessels of the lung pulmonary arteries and airways, evaluated by assessing 
the multifractal spectrum inherent in the power-law behavior, wherein the spectra of non-fractal, 
mononfractal and multifractal branching processes can be assessed(23). 

Patterns of complexity in a Self-Organized Critical State. 

Structural complexity in a self-organized critical state is distinguished by the absence of a 
characteristic scale(ll,12) representing a design objective that minimizes entropy and energy 
dissipation for a given metastable state (13). In transport networks in nature, structural complexity 
translates into asymmetric branching patterns modeled via a two-scale Cantorian branching 
process (12,24). Disorder is patterned through symmetry-breaking properties, leaving a branching 
structure devoid of a characteristic diameter scale, a characteristic bifurcation design and a 
characteristic fractal dimension (12, 17-19). Order is reflected via the respective fractal scaling of 
the symmetry breaking properties in the branching hierarchy. For biological networks we first 
characterize fractal scaling of network asymmetry by a self-similar fractal dimension based on 
bifurcation design x(3). A self-similar multifractal spatial distribution forms the basis for 
breaking the symmetry of bifurcation design within a tree(17-19, 23). However, in contrast to 
conventional fractal design theory, where scaling principles are assumed to be strictly self-similar, 
fractal design in fluid transport networks residing in a self-organized critical state include self- 
affine scaling(12, 14). With self- affine scaling, the fractal and multifractal dimensions change 
with branching generation when viewed from different perspectives, such as global and local 
regional scales(14, 16). This latter condition is of consequence to biological networks. At fine- 
grained diameter scales, changes in slope imply regional adaptations or mal- adaptations in 
branching design, not evident when viewed at a coarse-grained diameter scale. While these 
conditions formalize complexity in a branching structure in a critical state, and lead to 
complicated power-law behavior, they are sufficiently general to include self-similar fractal and 
non-fractal theories of universal design as special cases. 

Bifurcation design. 

Transport vessels carrying a flow, Q, branch by asymmetric bifurcation and possess a 
local symmetry-breaking architecture with a design related to their diameter. In a bifurcation the 
design, x, of the vessels is related to branching architecture and geometry on the basis of a 
relationship between parent (d 0 ) and daughter diameters (d h d 2 ) where , d 0 > >d 2 , 

< = < + d x 2 (1) 
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In biological vessels, such as arteries, diameters adhere dynamically to the condition Q oc d x 
where flow conservation condition Q 0 = Q ] + Q , holds (3, 20). Bifurcation architecture is 
conventionally summarized in terms of a diameter ratio, R d = d 0 /d { and asymmetry ratio, 
y = d 2 /d { . Scaling exponents, x, measured in networks characterize alternative theoretical design 
principles for a bifurcation which govern its fluid dynamic efficiency(3). The possible range of x 
is 0 < x < oo . However, the measured physiological range is smaller and include values of x=2, 3, 
and 4(1- 3, 20). In theory, these values form distinctive geometrical constraints for conserving or 
optimizing geometry and/or fluid dynamic properties in parent and daughter vessels (3). Design 
conservation principles between parent and daughter vessels occur at certain values of x, with x=2 
matching surface area, x=3 represents Murray’s law (1), here matching fluid volume; and x=4 
matches fluid conductance.. Under conditions of laminar flow these designs have potential 
physiological significance in that they conserve fluid velocity (x=2), shear stress plus hydraulic 
power (x=3) and pressure drop (x=4). In mammalian arteries, shear stress is considered an 
essential control element of design in both fractal and non-fractal theory(l, 3, 20) where the local 
coupling between shear stress as a threshold-dependent force and wall-stress as a threshold- 
independent force may constitute elements of two feedback mechanisms leading to a critical state 
(12, 20). The condition x=3 represents a stable flow condition(l), where energy dissipation and 
shear stress with flow are minimum. For x<3, flow is marginally unstable, with potential 
amplification of shear rates in daughter vessels relative to the parent(24). For x<2, an unstable 
flow condition, the area ratio of the bifurcation is less than one, and both flow velocity and shear 
stress are amplified in the daughter vessels. In this schema, maintaining constant shear stress 
requires the limitation of flow to below a critical threshold (9). 

Fractal dimension based on bifurcation design 

Bifurcation design has been thought as forming the basis of a fractal dimension in 
branching trees as generations increase to infinity(7). Asymmetric branching in biological 
networks (25) is analogous to a two-scale Cantorian branching process (23), with the parent 
diameter, d 0 , of equation 1 as the initiator and the two dissimilar daughter diameters represent the 
generators of the process (Figure 1). In generation n, there are N=2 n vessels. Branching to the left 

results is the major diameter pathway where vessel diameters scale according to R d n , while the 
minor diameter branching on the right, scales with each generation according to gR d ' n . Within a 
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generation there are 




vessels with diameter R d k (yR d l ) n * spanning from k=0,l,...,n. To 


formally define the fractal dimension, a measure M x is constructed, where 


M 


k= 0' 


O -foe 
\l 


(rtf * -(*? +(«■)' 


( 2 ) 


As n — > oo the diameter scale (yft , 1 j — > 0, such that M x remains finite if and only if x=D () , where 


D () satisfies 

«;°-+(* i r = i o) 

This limit defines a network design D 0 based on a bifurcation design x in terms of a Hausdorff- 
Besicovitch dimension. This branching process is said to be monofractal, lacking a characteristic 
diameter with generation, but exhibits order with branching through the fractal scaling 
relationship, D () , R d and yR ~ l . While x is an index conveying a design related to fluid dynamic 

efficiency (3), its convergence to a fractal dimension should not be confused with or identified 
with a topological dimension related to space-filling properties, except in the limit when 
branching design becomes space-filling(7, 15-16). 

Asymmetric fractal branching leads to a power-law relationship between a branching ratio 
R b , and the diameter ratio R d (25) The number of branches, N, as a function of diameter d, with 

each generation, g, is a self-similar recursive function 

N,(d) = N,_, ( dR ) • (4) 

D, introduced by Horsfield (25), accounts for the degree of asymmetry in numbers between 
daughter vessels and takes on integer values, such that D=0 reduces to a symmetrical branching 
tree and integers D>0 introduce successively greater degrees of asymmetry with branching. A 


branching ratio, R h is determined from a single iteration of the branching process if the limiting 
diameter at the termination is known to be fractal. The ratio of successive values of N , is 
R b = N JN t and when substituted into Eq. 4, gives R h ' (1+A) + R b l = 1, with one real positive root 

for R b . Denoting, R M =R b and R m = R' m a , such that ( R m > 2 ) and ( 1 < R M < 2), we obtain two 


alternative forms of the same underlying power-law 


Rm = & 


p(l+A)Z) 0 


R.„ = R 


( 5 ) 

(6) 


where y = R tl A and Rj + R m ' =1. Equation 5 and 6 are alternative scaling interpretations of self- 


similarity in power-law behavior dependent on a mono-fractal bifurcation design, x=X=D 0 . The 
slopes of the power-laws are both equal to x=D 0 but indicate different diameter scaling along the 
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major and minor branching pathways. Equation 5 defines a fine-grain diameter scaling of the 
power-law, where R d leads to an R M such that (1 < R M < 2). Equation 6 defines a coarse-grain 
diameter scaling of the power-law, where yR d 1 > R d yields R m > 2 . These branching relationships 
account for linear self-similar power-law behavior only, whereas power-laws with log-periodic 
branching can be observed. Such alternative behavior can be accounted for by a complex-valued 
fractal dimension(27). 

Conversely, self-similar power-law behavior in biological networks corresponds to a 
fractal branching process having two network interpretations. Eqs. 5 and 6 correspond to two 
classes of equivalent branching networks (14), which imply different connective properties in 
shaping the forces involved in self-organization (4). Fine-grained diameter scaling, (R d with R M < 

2) represents the connective branching properties of an equivalent asymmetric bifurcating 
network with a constant D(27). Alternatively, coarse-grained diameter scaling (g R rl 1 with R m > 2) 

represents the branching properties of an equivalent symmetric non-bifurcating branching network 
without connectivity, possessing fewer generations than the related asymmetric network for the 
same range of diameters (14). There are several methods for summarizing the connective and 
non-connective branching properties, R d , R m , R M and X=D () in biological branching networks. 
Such methods are useful, but artificial, being based on different topological approaches for vessel 
classification, such as the Horsfield (R V1 ~ R n < 2) and the Strahler (R m ~ R s > 2), or on ranking 
schemas(26). If the underlying branching process is a strictly self-similar branching fractal, then it 
does not matter what method is used to characterize X=D 0 (14). However, biological branching 
networks do not branch indefinitely (7), nor necessarily with the same bifurcation design filling 
space(4, 7). Therefore, the nature of the power-law behavior is not sufficient evidence of a fractal 
condition(15, 16). Also necessary for establishing a fractal condition is either evidence of self- 
similarity in scale-invariant linear power-law or a self-affine scaling condition in a non-linear 
scale-varying power-law(14). Such conditions are evaluated via a convergence test, accomplished 
by the evaluation of the multifractal spectrum (Figure 1) which tests for a fractal or non-fractal 
condition in the limit as branching becomes infinite (23). 

M ultifractal basis of bifurcation design. 

A mono-fractal network is not necessarily an optimal form that minimizes energy 
dissipation in fluid transport 17- 19) and does not account for the observation that bifurcation 
design in biological networks exhibits a distribution of values(4-5, 8). An additional symmetry- 
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breaking form of structural organization found in branching transport networks in non-living 
systems in a self-organized critical state includes mutifractal scaling(17-19). The branching 
structure of Eq. 2 becomes multifractal under slightly different scaling conditions, leading to a 
tree lacking a characteristic bifurcation design within a generation, but manifests order via a 
spatial distribution of bifurcation design exponents, analogous to x, defined in terms of a 

Hausdorff-Besicovitch dimension(23). Symmetry-breaking of the bifurcation design in the 
asymmetric mono-fractal network is introduced by an additional scaling relationship between g 
and R d , apart from the Horsfield scaling condition (25) R m = R' m a such that y = R d A , where A > 0, 

is a constant, here called a Woldenberg delta(28), where A A. A multifractal distribution is 
classically expressed in terms of a distribution of fractal dimensions with statistical moment q, 
ranging from -» < q < oo in three related senses(24): as a distribution of critical exponents t(q), 

as a generalized set of fractal dimensions, D q ; and of interest here, a multifractal spectrum, f(a) vs 
a, which relates to the spatial distribution of a, where a, is analogous to singular spatially 
determined values x, and f(a) corresponds to its density distribution. 

In a multifractal branching network bifurcations posses critical exponents x=t(q) where 

pX {q) + pX {q) = l (7) 

such that t(q) is a member of an underlying distribution function indexed by the moment q, Here, 
1, and 1 2 (l 2 >li) are the diameter contraction ratios ( /, = yR d 1 and l 2 = R~ d 1 ) in the minor and major 
daughter branches of a bifurcation where R d and gare defined above. The terms p , , p 2 (p 2 >pj) are 
defined in terms of R m and R M , where p i = R~ l <\, p 2 = 1 - p x = R~^ > \ , and are associated with 
diameter contraction ratios 1, and 1 2 . If Equation 7 is considered a generator in a branching 
process, we can subsequently define a measure in a manner similar to Eq. 2 above, 

m, = (pm + />,%')” = gf" = i (8) 

It is clear that as the number of generations n — > °° , M x remains bounded for the choice of 
x=t(q), M x — > 0 for a bifurcation scaling exponent x > r(c/), while M x —> for v < z(q) . This 

condition makes t(q) a critical exponent. In the multifractal formalism, t(q) is related to a 
distribution of fractal dimensions D = r(<y) /( 1 - q) for an infinite set of moments q over 
-oo < q < oo, bounded by D x > D 0 > D r The distribution D q has characteristic fractal 
dimensions related to specific moments q: for q=0, D 0 is the fractal dimension of design support 

as in Eq. 3; for q=l, D, is the information dimension; and for q=2, D 2 is the correlation dimension. 
Our interest is in another distribution called the singularity spectrum, f(a), which is the Legendre 
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transform of t(q), that is f(a) = min^a + r(<?)) and r(g) = max(/(a) - aa) . The value of q 
which minimizes qa+t(q) is the one in which a = -t'(<z) where the prime indicates the 
derivative. Here, q serves as a natural parameter indexing a and f in which a is bounded 
between a min = l) and a max = D_ oo . Now, the parameter a is an exponent measuring the strength 

of singularity or degree of unbalance of branching at a location within the branching process. The 
expression f(a) is the corresponding fractal density, indicating how large a part of the whole 
branching tree possesses this singularity of degree a, where the singularity strength of a within a 

generation n is approximately 2" ,(</| . 

Figure 1 illustrates multifractal spectra for a monofractal branching process compared to 
multifractal branching process. The multifractal branching network pattern illustrates one possible 
configuration with a simple form of spatial clustering of scaling exponents a(e) within the major 

and minor pathways leading to the indicated convex spectra. The spatial properties of a are 
specified as a position index e , wheree = and 0 < e < 1 in the branching process of Eq. 8. 
Here, a(e) are diameter scaling exponents analogous to x in Eq. 1, but are now singular values in 
a multifractal distribution whose values range from a min to a max , and whose value within the 
network depends upon bifurcation position defined by ein the binomial expansion of Eq. 8. The 
density spectrum of a ( e Is a limit probability distribution, f(e), as the order of branching 


becomes large, 


_ m 

(1 + Ae)ln R d 


(9) 


such that H(e ) = elne + (1 - e)ln(l - e) is the entropy distribution function of statistical 
mechanics and is convex in shape reaching a maximum of 1 at e=l/2. The relative magnitude of 


A and D influences the range of extreme values of a 


a( 0)= lni?M 
lnfl/ 

(10) 

a( 1) = l + _a(0) 
v 7 1 + A v 7 

(11) 


The density f(a(| )), corresponding to e = l/2i§ simply the fractal dimension of the support related 


to the moment q=0, specified by f(a( \ ))=D 0 . 


A power law does not convey whether an underlying branching process is multifractal and 
the power-laws implied by a branching process converging to a convex multifractal spectrum are 
not simple. The generalized fine-grained and coarse-grained power laws in multifractal branching 
are of the form 
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(l+Ae)/(l+Ae) _ a (£) 

~ n d 


_ a(s)( 1+Ae) 

~ 'Vi 


( 12 ) 

(13) 


where symmetry-breaking emerges under conditions of A ^ A but reduces to a mono-fractal 
branching architecture when D= A , and where D and A no longer are restricted to integer values. 

There are an infinite number of power-laws according to the range of e, but power laws analogous 
to Eq. 5 and 6 are obtained by setting e=l/ 2where the slope of the power-law is equal to D 0 , the 


same result for the monofractal branching condition. However, unlike monofractal branching, a 
multifractal configuration imparts a spatial distribution of bifurcation cross-sectional area capable 
of modulating the fluid dynamic forces involved in self-organization(29). The power-laws in Eq. 
12 and 13 are linear but can have log-periodic variations if D 0 is a complex-valued fractal 
dimension(27), or be self-affine(14). Consequently, a multifractal branching condition is verified 
via the computation of the multifractal spectrum and its convergence to a convex spectrum(23, 

30). 

Self-affine fractal scaling of bifurcation design 

An additional property serving to break symmetry in a self-organized critical branching 
network is to modify the form of fractal scaling in the branching process, from a self-similar 
scale-invariant state, as assumed in present simple-fractal branching theory, to a self-affine form 
of scaling, which is scale-varying(12, 14). Under these circumstances, the fractal dimension D 0 
varies continuously with the branching region of interest. Consequently, non-linear power-laws 
arise, one form which may appear piece- wise linear, under fine-grained diameter scales, a 
property found in the power-laws of lung airways and pulmonary arteries(14). A verification of a 
fractal self-affine condition, as opposed to the branching process not being fractal, is established 
by evidence of dissimilar convergent multifractal spectra over different regions of branching 
where the slope of the local power-law becomes the local regional fractal design of support, 
D 0 (23). 

A departure from a self- similar scaling condition leads to locally emergent behavior in 
branching properties for biological networks. Self-affine scaling of bifurcation design is expected 
and necessary for a physiological branching network to adjust its bifurcation design, x, toward a 
value of 3 as branching becomes space-filling in exchange vessels (7, 16). In a self-affine 
condition, the global scaling properties of fractal design, which characterize the energy dissipation 
in a marginal stable state (12) are now disassociated from the scaling properties responsible for 
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the modulating the distribution of forces involved in self-organization(4). This disassociation 
implies that spatially varying designs are present throughout a branching tree while the global 
power-law appears unchanged at a coarse-grained diameter scale. Consequently, self-affine 
scaling potentially identifies hidden regional network scaling adaptations at fine-grained diameter 
scales, where fractal design shifts away from ideal bifurcation behavior of x=3 or ideal power-law 
behavior of X=3. Also, self-affine scaling eliminates the possibility that a branching network 
exhibits a universal design with characteristic energy dissipation properties(l,7). 

Complexity versus Universal design in the lung 

While one cannot mathematically prove fractal behavior in limited branching networks 
(15), one can alternatively assess conclusive evidence of the presence or absence of a 
convergence criterion as branching becomes infinite via the evaluation of the multifractal 
spectrum(23, 30). Therefore, testing for fractal complexity versus universality was evaluated by 
comparing power-laws and multifractal spectra in the lung airways (31) and developing 
pulmonary arterial system (8,32) using fine-grained and coarse-grained diameter scales. 
Complexity in the fetal/neonatal pulmonary arterial circulation 

One of the most dramatic physiological adaptations made by a biological transport 
network in nature occurs at birth in the mammalian pulmonary vascular system after the first 
breath (10, 33). The lung in the fetus is a dormant organ with a high resistance, requiring a small 
level of blood flow for development. Arterial vessels have made structural adaptations to a 
physiologically high pressure that is never again experienced postnatally, except under extreme 
exercise conditions, or disease states(34). After the first breath of life, lung airway expansion is 
accompanied by an arterial network transition from a high-pressure, low-flow network to a lower- 
pressure, high-flow system( 10, 33). While the distribution of x in the fetal and neonatal 
pulmonary circulation is unknown, Thompson argues that a starting working hypothesis would be 
based on the principle of minimum work in bifurcation design, where x=3(l, 9). Consequently, 
this would imply that the high resistance to blood flow in the fetal arterial system is dictated by 
humoral factors constricting vessel diameters, whereby the diameters increase with flow and 
shear-stress dependent factors along with lung expansion after birth(33). However, a phenomenon 
not explained at this stage of development, is an unusual susceptibility to flow-induced vascular 
injury when relatively modest elevations in blood flow above a threshold occur, yielding 
additional structural modifications in arteries leading to complications in adaptation at birth(34). 
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Hopkins theorized that structural or mechanical factors within the pulmonary circulation cause the 
perinatal arterial network to posses a different profile of energy dissipation than that of a mature 
pulmonary arterial system(35). A corollary hypothesis is a self-organized critical interpretation, 
predicting that the arterial design is in a metastable state, possibly far removed from a global 
minimum of energy dissipation of x<3. 

Figure 2 represents three theoretical views of bifurcation design in a fetal lamb (140 
days/145 days at term) Fig 2ACE, and a 1 day-old neonatal lamb. Fig 2BDF, representing the 
structural state, just prior to, and just after birth(8). Figure 2A and 2B represents the distribution 
of bifurcation design exponents, x, for the fetus and newborn, calculated for each bifurcation in 
the lung calculated according to Eq. 1. The gray area in Fig 2AB and 2EF marks those 

bifurcations with a design condition of x<2, an unstable condition for flow transport, indicative of 
an area-ratio, b<l. The figure indicates that a significant proportion does not adhere to ideal 

principles, but does not reveal where they are located or how they are scaled with branching. Fig 
2CD represents the corresponding power-laws at a coarse-grained (squares) and fine-grained 
(circles) scaling of diameter. The red, blue and green colors correspond to the slopes of regional 
power-laws indicative of the condition x=X in different regions of the lung relative to the main 
pulmonary artery feeding the arterial system: globally, proximally and distally, respectively. The 
changes in slope do not identify whether the system is non-fractal, or possesses a complex-valued 
fractal dimension (28), or if the fractal branching is self-similar or self-affine. The fine-grained 
diameter scaling power-law in Figure 2CD suggests that the location of bifurcations with values 
of x<2 in Fig 2A and 2B are located distally as one branches outward towards smaller vessels. 
Figure 2EF represents the evaluation of the fetal and neonatal multifractal spectra and the 
emergence of complexity from Figure 2AB and 2CD. The colors in Fig 2EF mark regions of 
branching corresponding to the power-laws of C and D where the slopes of the power-laws are 
equal to the peaks of the corresponding multifractal spectra. As a convex spectrum represents a 
multfractral branching process as one approaches an infinite branching limit, one concludes that 
the underlying distribution of bifurcation design is multifractal, as such a condition would not 
have arisen had the branching process been simply stochastic or strictly self-similar being 
characterized by a unique dimension(23). That a different multifractal condition is present in other 
branching regions implies that the branching process is self-affine, not strictly self- similar! 14). 


12 



These results indicate that in the higher pulmonary resistance state of the fetal and 
neonatal pulmonary circulation, the larger pulmonary arteries do not follow universal biological 
design principles of x=2 or 3, instead exhibit fractal branching complexity in a metastable state of 
marginal stability far removed from a global minimum of energy dissipation(12). After birth, 
there is an adaptation of design, regions shifting globally and locally toward higher values, but 
remain far removed from a universal minimum of x=3. It is important to emphasize that at a fine 
grained level of diameter scaling, the proximal and distal regions in the fetus and newborn show 
disparities in design with the distal arteries with a more severe decrease in cross-sectional area. 
The design in the distal region of both fetus and newborn represents a bottleneck of decreasing 
cross-sectional area with branching, where shear stress associated with blood flow is amplified as 
it progresses toward the periphery(4). Above a threshold, shear stress associated with blood flow 
is capable of inducing a chain of events occurring over different time-scales, leading to changes in 
the force-transduction properties of the endothelial wall(20) and the integrity of vessel structure 
and function, with large increases in flow consequently producing injury(34). The unusual 
bifurcation design of x<2, along with the branching complexity inherent in the structural 
organization of the fetal and neonatal pulmonary circulation therefore explains a basis for 
Hopkin’s theory of dissipation and for the unusual susceptibility of the perinatal pulmonary 
circulation to injury from elevated blood flow(35) 

Complexity of Adult Human Pulmonary Arterial Circulation 

While the developmental patterns of arterial branching complexity in the pulmonary 
circulation are not known, adulthood is expected to adhere to a universal design (1, 2, 7) as an 
adaptation of design that is more efficient than the fetal and newborn state (10). However, while a 
universal design prescribes an allometric scaling of resistance at a fixed flow rate commensurate 
with body weight (7), such a theory is inadequate in predicting the variations in resistance 
between individuals over a wide range of flows, as experienced at rest and extreme exercise(36). 
For an unknown reason, some individuals, including highly trained athletes, have highly elevated 
resistances above the expected universal norm and are fatally susceptible to consequences of 
pulmonary hypertension under extreme environmental and exercise conditions (36-37). 

Analogous to the fetus and newborn pulmonary circulation, the susceptibility is hypothesized to 
have structural and/or humoral origins(35). A universal design hypothesis for the pulmonary 
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circulation implies that a self-similar fractal dimension will be observed (7) while a complexity 
hypothesis can implicate a design variation that is apparent only at fine-grain diameter scales. 

Figure 3ACE represents a fractal analysis of an adult human pulmonary vascular system 
free of pulmonary vascular disease (32) expressed in the same format as Figure 2. In Fig 3A, the 
distribution of bifurcation design exhibits a large proportion of values less than 2, consistent with 
local decreases in bifurcation cross-sectional area with branching. Fig 3C demonstrates that, 
similar to the fetus and newborn state, regional self-affine scaling of fractal design is evident 
when viewed from fine-grained diameter scales. In Figure 3E, the convergent multifractal spectra 
with different fractal dimensions of support are found in different regions of arterial branching, 
indicating that the network branching design scales in a self-affine manner. In Fig3E, the distal 
branching indicates a regional instability of decreasing cross-sectional area where x<2, a 
condition similar to that found in the fetal and newborn state. Such a decrease in cross-sectional 
area, not detectable from a coarse-grained view of diameter scaling represents a persistence of the 
fetal design and a possible structural origin of susceptibility, analogous to the fetal condition of 
blood flow injury secondary to elevated shear stress and blood flow(354-35). 

Complexity of Adult Human airways 

The laws of lung development (38) dictate that the airway branching pattern is laid down 
prior to the arterial vascular branching pattern, and the latter arises from the airways before blood 
flow and shear stress can contribute to such a process (7). However, as the previous results 
suggest, arterial design adapts dramatically during development and growth, but airway branching 
power-laws are consistent with Murray’s law of x=X ~ 3 at birth (39) indicating a different 
process of design development. Under a universal design hypothesis, a self-similar monofractal 
design of x=3 expected to be present. Alternatively, under a self-organized critical hypothesis, the 
complexity of the airways is expected to give rise to complexity in the arterial system where both 
lack a characteristic fractal dimension of design. 

Figure 3BDF summarizes the branching analysis for the airways. Fig. 3B illustrates a 
greater proportion of bifurcations around x=3. Figure 3D shows linear power-law behavior close 
to x=X=3 for coarse-grained diameter scales, but with additional slope changes with fine-grained 
scaling. Figure 3F, is the multifractal analysis of airway branching data of an adult human. This 
figure demonstrates that the airways exhibit branching complexity, with multifractal branching 
and a self-affine state of fractal scaling, consistent with a system in self-organized critical 
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state(12). While the coarse-grained diameter scaling indicates a design consistent with Murray’s 
law, the fine grained diameter scaling indicates that smaller vessels change their fractal dimension 
of support with branching. This result shows that the complexity of bifurcation design, laid down 
early in development, is apparent in the adult airways, and may therefore influence the 
development of design in arterial pathways as well(10, 38). 

Conclusions 

The design of larger transport vessels in the lung were found to exhibit fractal branching 
complexity with an origin consistent with systems in nature conforming to a self-organized 
critical state. The changes complexty in the pulmonary circulation are dynamic suggesting that 
the lung is a complex adaptive system lacking a lacking a universal fractal design, and is not a 
strictly self-similar branching network. Also, there is evidence of a fractal self-affine form of 
scaling, implicating a structural organizational source of vascular disease susceptibility. However, 
while such complexity appears to account for much of the bifurcation design variation in the lung, 
it is not apparently responsible for all possible patterns seen in branching morphogensis (40), 
therefore, should not be considered the only process governing biological branching pattern 
formation. Alternatively, a self-organized critical state may be part of a universal pattern of 
biological network development, growth and adaptation. Even with this limitation, fractal 
branching complexity of a self-organized critical network may obviate the concept of a universal 
design (1, 2, 7), remaining consistent with universal biological design principles (41) while 
providing a theoretical framework for understanding the process of biological network adaptation 
(9,10). 
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Figure 1. Architecture of a branching network. Alternative theories of biological network design predict different 
consequences for the scaling and distribution of bifurcation design x in relationship J 0 ' = d\ + d* in trees. 

Consequences manifest themselves in the spatial distribution of x within branching tree, and in the appearance of 
multifractal spectrum, f(a) versus a, which tests convergence to a limit probability distribution of branching structures 
with limited generations as branching becomes infinite. Non-fractal theories predict x is distributed randomly without 
scaling relationships, monofractal theory predicts identical distribution of x in tree. Self-organized critical state predicts a 
multifractal scaling relationship spatially distributing x in network along with a self-affine scaling of the fractal dimension 
with generation. 

Bifurcation design x: Between parent and daughter vessels, the exponent x dictates fluid dynamic efficiency of fluid 
transport and cross-sectional area with branching: x=2 is constant fluid velocity with the preservation of cross-sectional 
area; x=3 is Murray’s law corresponding to constant fluid volume with constant hydraulic power; and x=4 represents 
constant fluid resistance with preservation of pressure drop. For values x<2 cross-sectional area of bifurcation decreases, a 
condition resulting in amplification of flow velocity and shear stress, forces participating in vessel self-organization. 

Multifractal branching condition represents one aspect of branching complexity inherent in a self-organized critical state. 
Branching condition involves symmetry-breaking of mono-fractal pattern into a spatial distribution of scaling exponents 
such that X — > a(e), where a (e) is a singular distribution of bifurcation exponents distributed with branching according 
to address e, where 0<e<l, Density of a ( e ) is given by £( ( e ) ) in multifractal spectrum. Conditional) marks bifurcation 
exponent a(0) along major daughter pathway e=0, and e=l marks bifurcation exponent a(l) along minor daughter 
pathway. Fractal dimension of support for bifurcation design is designated by £ = j where Oi(^) = D () . Address of 
bifurcation a ( e ) in tree is determined by index k and n of binomial expansion of Equation 7 in text. In figura ( e ) is 
distributed over a range of values a min to a max . encoded with colors of rainbow. If a min =a(0) of major diameter pathway, 
then minor pathway must possess condition a max =a(l) and vice-versa. Control of distribution of scaling exponents a (e) and 
shape of multifractal spectrum, is via the terms A and A , called the Horsfield and Woldenberg delta in text. 

Self-similar versus Self-affine fractal scaling of bifurcation design: The scaling of D 0 with branching is different in the 
three theories. Non-fractal theories do not exhibit scaling properties, hence no fractal dimension. Universal fractal theories 
presuppose the fractal scaling is self-similar, leading to conditon of scale-invariance of design throughout branching tree 
and strictly linear power-laws. Complexity interpretation predicts multifractal branching plus self-affine fractal scaling 
whereby the D„ varies from one generation to another leading to non-linear power-laws. 
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Figure 2. Fetal and Neonatal Ovine Arterial System. Three scaling views of bifurcation design in fetal and neonatal 
pulmonary circulation where branching architecture does not adhere to idealized design principles of x=3: AB -distribution 
of x as if a stochastic branching network. CD-power law behavior for different regions of branching, EF-multifractal 
spectra of branching as a network in a self-organized critical state. Data from (8). 

AB -Distribution of bifurcation design in tree: Bifurcation design, x, in each bifurcation solved according to 
cIq = d'l + d] where 0 is parent diameter, 1 is major daughter and 2 is minor daughter. Area ratio, b is sum of squares of 
daughter vessels divided by square of parent area. Grey area marks bifurcations with area ratio<l. Fetus (N=2844 
bifurcations where bifurcation = 3 vessels). Fetal mean values of x: arithmetic mean=2.076; geometric mean=1.684; 
harmonic mean =1.50. Fetal mean values b: arithmetic mean=0.8573; geometric mean=0.81 1; harmonic mean 0.753. 
Newborn (N=812 bifurcations). Newborn mean values of x: arithmetic mean=3.44;geometric mean=2.059; harmonic 
mean=1.453. Newborn mean values of b: arithmetic mean=0.944; geometric mean=0.767; harmonic mean 0.469 

CD- Power-law behavior of distribution in AB: Underlying distribution of bifurcation exponents exhibits a non-linear 
power-law behavior suggesting that branching is a self-similar fractal with a complex-valued fractal dimension predicting 
log-periodic behavior that is scale-invariant. Non-linear interpretation draws slopes in proximal and distal regions for sake 
of comparing with global design. Coarse-grained (R dc =squares) and fine-grained (R d[ =circles) decimation of diameter in 
bifurcations were used to form corresponding power-law relationships, R" n = R"* and R" M = R'^ with equal slopes = X 

(colored red) where n denotes generation starting from n=0 for largest vessels inceasing in index towards smaller vessels. 
Coarse graining with diameter scale R dc yields an equivalent symmetric branching network without bifurcation connectivity 
having a branching ratio R m =R b >2. Fine-graining with diameter scale R df yields equivalent asymmetric branching network 
with bifurcation connectivity having R M =R b <2. Fetus: branching properties R dt =1.33, :R dc =2. 1 52; R M =1.51; R m =3.0. 
Newborn: branching properties R df =1.32, R dc =1.613, R M =1.66, R m =2.52. Fetus: values of power-law slopes X: Red- 
global X = 1 . 4 3 ; ft /«c - p ra x i m a 1 = 1 . 9 94 ; G ree n - d i s t a I =0 . 9 3 9 . Newborn: values of power-law slopes X: Red-global 
= 1 ,79;B/Me-proximal=2.44;Green-distal= 1 .40. 

EF -Multifractal Spectra of power-law behavior in CD: Convergence to convex multifractal spectrum in proximal and distal 
regions evidence of self-affine fractal scaling with each local D 0 equal to local slopes X. A global view of branching is 
colored red with superimposed multifractal spectra: symmetric branching-coarse grained(boxes) and asymmetric -fine 
grained (circles). Proximal region colored blue and Distal region colored green. While regions were chosen arbitrarily to 
compare larger versus smaller vessels in branching, any combination of branching regions chosen manifest multifractal 
branching designs, indicating that the fractal dimension of support D 0 scales continuously with generation. Results support 
theory that branching architecture is in a self-organized critical state where complexity in branching yields a global design 
minimized for energy dissipation but varies local design without a characteristic design, to modulate forces of self- 
organization throughout tree. 
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Figure 2. Ovine Fetal and Neonatal Pulmonary Arterial Systems 
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Figure 3. Human Pulmonary Arterial (ACE) and Airway System (BDF). Pulmonary arterial system branching is 
known to arise from branching morphogenesis of airways. Stochastic bifurcation design theory and mono-fractal theory 
predict fixed designs at x=3 or x=2 in larger transport vessels. Theory of Self-organized critical state predicts complexity of 
arterial system is found in complexity of airways, but not necessarily with the same design or scaling. 

Arterial System -(ACE)-data from S.S Singhal, R, Henderson, K. Harding and G, Cumming Circ Res 33 190 (1973). 

N=893 bifurcations 

A -Distribution ofx and b . Grey area marks bifurcations with area ratio b<l. Mean values of x: Arithmetic mean=2.33, 
geometric mean =2.23; harmonic mean=2.12. Mean values of b; arithmetic mean=1.063; geometric mean=1.045; harmonic 
me an= 1.025 

C -Power-law behavior of distribution in A: As in Figure 2, coarse-grained (R dl .=squares) and fine-grained (R d ,=circles) 
decimation of diameters in bifurcations were used to form corresponding global power-law relationships, R" t = R"* and 
R" m = Rjj with equal slopes = X (colored red)=2.39. Global arterial branching properties R df =1.33, :R dc =2.152; R M =1.51; 

R m =3.0. At fine-grained decimation, power-law exhibits changes in power-law slope suggestive of a change in fractal 
design. Value of slopes; X=Proximal branching-(blue)=2.62, X=Distal (green)=1.81. Distal region implicates decreasing 
cross-sectional area with branching where non-ideal bifurcation designs reside as a population, similar to what was 
observed in fetus and newborn lamb. 

E -Multifractal Spectra of power-law behavior in C: Convergence to convex pectrum globally, proximally and distally, is 
evidence that underlying fractal dimension is self-affine multifractal consistent with structure in a self-organized critical 
state. The global view of branching (red) as a multifractal network can also be viewed as two branching processes with 
successive fractal designs proximal and distal. The distal design does not adhere to ideal design principles and appears to 
be a developmental remnant of fetal state. 

Airway System-{ B DF)-data from K. Horsfield, J Appl Physiol 24 373 (1968) N=2613 bifurcations 

B- Distribution ofx and b . Mean values of x; Arithmetic mean=2.78, geometric mean =2.53; harmonic mean= 2.35. Mean 
values of b; arithmetic mean=1.124; geometric mean=1.102; harmonic mean=1.079 

D -Power-law behavior of distribution in A: As in Figure 2, coarse-grained (R dl .=squares) and fine-grained (R d( =circles) 
decimation of diameters in bifurcations were used to form corresponding global power-law relationships, R m = R* c and 
R M ~ R d f with equal slopes = X (colored red)=3.1. Global arterial branching properties R dl =1.157, :R dc =1.387; R M =1.59; 

R m =2.7. At fine-grained decimation, power-law exhibits changes in power-law slope suggestive of a change in fractal 
design. Value of slopes; X=Proximal branching-(blue)=3.0, X=Distal (green)=2.51. Distal region has significantly different 
design of power-law behavior than proximal region. 

F -Multifractal Spectra of power-law behavior in C: Convergence to convex spectra globally, proximally and distally, in 
airways is evidence that underlying branching network is a self-affine multifractal consistent with structure in a self- 
organized critical state. The complexity of the airways appears related to the complexity of the arterial branching. 
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